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Abstract 

The mechanism of color confinement as a consequence of an unbroken non-abelian 
gauge symmetry and asymptotic freedom is elucidated and compared with that of 
other models based on an analogy with the type II superconductor. It is demon- 
strated that a sufficient condition for color confinement is given by Z^ 1 — where 
Z3 denotes the renormalization constant of the color gauge field. It is shown that 
this condition is actually satisfied in quantum chromodynamics and that some of 
the characteristic features of other models follow from it. 



1 Introduction 

It is our consensus that strong interactions are governed by quantum chromodynamics 
(QCD) or the gauge theory of quarks and gluons. These fundamental constituents of 
hadrons carry the color quantum number and are considered to be unobservable. This 
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is a conclusion drawn from our unsuccessful attempts to observe isolated quarks and is 
referred to as color confinement - abbreviated as c.c. - in what follows. It is the subject 
of this paper to compare various interpretations of c.c. in an effort to extract common 
features from them. For this purpose we review early attempts to interpret confinement 
in Section 2. These theories are formulated in configuration space and are characterized 
by the two phase structure or the dual Meissner effect with a finite penetration depth. In 
Section 3 we recapitulate an interpretation of confinement formulated in the state vector 
space on the basis of BRS invariance and asymptotic freedom (AF) skipping details of the 
proof. In this section we present two alternative conditions for confinement. In Section 4 
we prove that one of these two conditions is a consequence of the other on the basis of the 
renormalization group (RG) method. It is shown that c.c. is realized when the condition 
Z^ 1 = is satisfied where Z3 is the renormalization constant of the color gauge field. It 
is called the condition for color confinement abbreviated as CCC hereafter. In Section 5 
it is proved that the CCC is really satisfied in QCD. 

Finally in Section 6 we compare the consequences of the CCC with other interpreta- 
tions. First, we give an intuitive interpretation of the CCC in fictitious electrodynamics. 
Next, we discuss the connection between the CCC and the linear potential between a 
heavy quark and a heavy antiquark pair resulting from Wilson's area law. Then, we show 
that the flux of the color gauge field emerging from color singlet hadrons cannot penetrate 
into the confining vacuum leaving no trace of long range forces. This resembles the dual 
Meissner effect introduced in some of the other interpretations. 

Two appendices are included. Appendix A is intended to clarify the relationship 
between two alternative forms of the conditions for color confinement. Appendix B gives 
a derivation of (3.48) in the absence of asymptotic fields due to infrared singularities. 

2 Early Attempts to Interpret Confinement 

Because of its profound mysterious nature exhibited in strong interactions various attempts 
have been made to understand the mechanism of color confinement on the basis of classical 
or semi-classical gauge theories eventually exploiting topological quantization. 

Starting from a classical Lagrangian of the Higgs model, Nielsen and Olesen [1] identi- 
fied the hadronic strings with the Landau-Ginzburg-Abrikosov vortices of quantized mag- 
netic flux in the superconducting vacuum. Their vortices are either endless or closed and 
the energy of the system is minimized for a certain optimum radius of the vortex, and the 
total flux is topologically quantized. 

Nambu [2] introduced Dirac's monopoles into this theory and realized finite vortices 
by putting monopoles at both ends, and Dirac's quantization for monopoles matches the 
flux quantization. Thus a hadronic string is formed in a superconducting vacuum by 
joining a monopole-antimonopole pair by a vortex of a quantized magnetic flux. The flux 
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cannot spread out and the energy of the flux is proportional to the length of the string. 
This pair may be regarded as a meson and monopoles in the superconducting vacuum 
are then confined. An ordinary superconductor is a coherent superposition of charged 
objects, but a vacuum which confines electric, instead of magnetic, monopoles may be 
a coherent superposition of magnetic monopoles as suggested by Mandelstam [3] and by 
't Hooft [4]. Magnetic monopoles do not appear in continuum quantum electrodynamics 
(QED) unless they are put by hand, so that electric charges cannot be confined in QED 
if this interpretation should be taken for granted. They can appear in non-abelian gauge 
theories, however, as shown by 't Hooft [5], and therefore the Yang-Mills vacuum could 
be a coherent superposition of magnetic monopoles, and confinement must be a property 
characteristic of non-abelian gauge theories. We shall come back to this subject later again 
in Subsections 3.4 and 6.1. 

In the quantized version of non-abelian gauge theories, the only known way of including 
higher order corrections in a gauge-invariant manner is to exploit the lattice gauge theory 
[6], although Lorentz invariance is recovered only in the limit of the vanishing lattice 
constant. Since confinement is a non-perturbative effect it is important to check it in the 
lattice gauge theory. Wilson [6] has formulated confinement in the form of the area law for 
the loop correlation function leading to a confining linear potential between a quark and 
an antiquark. In the strong coupling approximation the area law is obeyed even in QED, 
but it is a non-trivial problem to check if the strong coupling regime could be continued to 
the weak coupling one without encountering a phase transition. Susskind and Kogut [7] 
analyzed the lattice gauge theory in the strong coupling approximation and found that the 
confining strong coupling phase resembles that of a type II superconductor with electric 
and magnetic fields interchanged. We shall come back to the discussion of the area law 
later again in the Subsection 6.2. 

All such attempts including the recent supersymmetric Seiberg-Witten theory [8] tend 
to indicate that confinement is a consequence of the coherent superposition of magnetic 
monopoles in the vacuum state dual to the superconducting one based on coherent su- 
perposition of charged objects such as Cooper pairs. Quarks and antiquarks are joined 
together by the electric flux penetrating into the vacuum or through the normal conducting 
phase enclosed by the superconducting phase. In other words, confinement is characterized 
by the two-phase structure or the Meissner effect of a finite penetration depth. 

In electrodynamics all the physical media are characterized by their dielectric constant 
e > 1. Suppose that there is a fictitious medium of dielectric constant [9], 

e<l. (2.1) 

This medium is antiscreening. When a small charge is placed in this medium, it will crack 
and develop a hole surrounding this charge. We have e = 1 inside and e << 1 outside. 
Because of the antiscreening nature of the medium the induced charge on the inner surface 
of the hole is of the same sign as the originally inserted charge. In order to reduce the 
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size of the hole the repulsion between the original charge and the induced one must be 
overcome, but elimination of the hole would require infinite energy, so that the hole will 
not disappear. The two media, one with e = 1 and the other with e << 1, are considered 
to correspond to normal conducting and superconducting states, respectively, and such a 
two-phase structure of the vacuum is common to all the models constructed in analogy 
with the type II superconductor in the configuration space. 

In the next section we shall recapitulate the arguments based on the BRS invariance 
[10] and asymptotic freedom [11,12] in the state vector space. 

3 Color Confinement as a Renormalization Effect 

In a covariant quantization of gauge fields introduction of indefinite metric is indispensable. 
Thus the resulting state vector space V involves unphysical states of indefinite metric and 
we have to find a criterion to select physical states out of V. For this purpose we employ 
the Lorentz condition in QED as a subsidiary condition, but it is more involved in non- 
abelian gauge theories. In what follows we shall confine ourselves to QCD, and in order 
to fix the notation we start from its Lagrangian density in the Pauli metric, 

£ = C-inv + C g f + £fp, (3-1) 

where 

Cinv = ~\F^u ■ F^ u - tp^^D^ + m)ip, (3.2a) 

Cgf = ■ d^B + %B-B, (3.26) 

C F p = id^c ■ D^c. (3.2c) 

We have suppressed the color and flavor indices above. The first term in (3.1) is the 
gauge-invariant term, the second one the gauge-fixing term and the last one the Faddeev- 
Popov ghost term. In (3.2b) a denotes the gauge parameter and B the Nakanishi-Lautrup 
auxiliary field, and in (3.2c) the hermitian scalar fields c and c are anticommuting and are 
called Faddeev-Popov (FP) ghost fields. 

In what follows we shall introduce the inner and outer products of two colored objects: 

S ■ T = S a T a , (3.3) 



(SxTT = Y,fabcS b T c , (3.4) 

b.c 
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where a, b, c etc. are color indices and f a b c the structure constant of the algebra su(3) 
corresponding to the color gauge group. Then covariant derivatives are defined by 



= (d„ - igT ■ Apty, (3.5) 

V = 9 /+^ xc )". ( 3 - 6 ) 
F- u = d^K - d u A- + g(A^ x A u ) a . (3.7) 

The local gauge invariance is respected by (3.2a) but not by the other two, (3.2b) 
and (3.2c), introduced for quantization. It so happens, however, that the total Lagrangian 
density is invariant under new global transformations called the Becchi-Rouet-Stora (BRS) 
transformations [10]. 



3.1 BRS transformations 

Let us consider an infinitesimal gauge transformation of the gauge and quark fields and 
replace the infinitesimal gauge function either by c or c. They define two kinds of BRS 
transformations denoted by 5 and 5, respectively. 



SA^ = D^c, SA^ = D^c, 
Sip = ig{c ■ T)ip, Sip = igic ■ T)ip. 



(3.8) 
(3.9) 



For the auxiliary fields B, c and c local gauge transformations are not even defined, but 
their BRS transformations can be introduced by requiring the invariance of the local 
Lagrangian density, namely, 



SC = SC 

then we find the following transformations: 



0, 



SB = 0,Sc = iB, Sc = —\g{c x c), 



SB = 0,Sc = iB,Sc 



C X c 



where B is defined by 



B + B - ig(c x c) = 0. 



(3.10) 

(3.11a) 
(3.116) 

(3.12) 



Noether's theorem states that the BRS invariance of the Lagrangian density leads to 
conservation laws and we introduce two conserved BRS charges Qb and Q B by 



i [Qb, ' 



(3.13) 
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where we choose the — (+) sign, when the field <fi is of an even (odd) power in the ghost 
fields c and c. 

Maxwell's equations for the gauge field can be expressed in terms of BRS transforma- 
tions as 

+ g.J v = i55A u , (3.14) 

where J u denotes the color current density and g the gauge coupling constant. This set 
of equations has been derived from the Lagrangian density (3.1), but the form of these 
equations is valid for a wider class of theories such as supersymmetric and grand unified 
theories as long as the original gauge symmetry is respected. The only deviation from 
the original theory represented by (3.1) appears when we try to express the color current 
density J v explicitly in terms of the elementary component fields. 

The r.h.s. of Eq. (3.14) represents a conserved current and stands for the deviation 
from Maxwell's classical equations. 

d^iSSA^ = 0. (3.15) 
The BRS charges are hermitian and nilpotent, for example, 

Qb = Qb, Ql = 0. (3.16) 

The nilpotency implies introduction of indefinite metric and a physical state |/) is defined 
by a constraint [14] 

Q B \f) = 0, |/)€V. (3.17) 

The collection of physical states including the vacuum state |0) forms the physical subspace 
of V denoted by V p h ys , 

V phys = {\f):Q B \f) = 0,\f)eV}. (3.18) 
It should be mentioned that the S matrix is BRS invariant, 

5S = i[Q B ,S] = 0, (3.19) 

so that the physical subspace V p h ys is an invariant subspace of the S matrix. 

Furthermore, we introduce a subspace of V called the daughter subspace Vd defined by 

V d = {\f):\f) = Q B \g),\g}eV}. (3.20) 
Because of the nilpotency of Q B , Vd is a subspace of V p h ys , 

V d C V phys , (3.21) 
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and we introduce the Hilbert space H by 

U = V phys /V d , (3.22) 
which may be called the BRS cohomology [15,16,17]. 

3.2 Relation to QED 

QED is the oldest example of gauge theories so that the above formulation should be 
applicable to it. We recognize that the subsidiary condition (3.17) to select physical states 
looks completely different from the Lorentz condition, 

£(+)(*)!/} = 0, (3.23) 

so that we shall clarify the relationship between them. 

In QED or in an abelian gauge theory the auxiliary fields B, c and c are free and 
massless, namely, 

UB(x) = Uc(x) = Uc(x) = 0. (3.24) 

Furthermore, in the conventional treatment of QED the ghost fields do not participate 
in the game so that we introduce constrained physical states [18] in terms of the positive 
frequency parts of ghost fields by 

Q B \f')=0, (3.17) 

c(+)(x)|/') = c(+)(x)|/')=0, (3.25) 

and the corresponding subspace V phys . Consistency among them requires the following 
condition: 

{QB,<t + Hx)}\f') = Bl+\x)\f')=0, (3.26) 

which is precisely the Lorentz condition. 

On the other hand, when we have the Lorentz condition and (3.25) the condition (3.17) 
follows automatically from the structure of Qb in QED, namely, 

Q B = Jd 3 x[B(-\x)dM{x) - B(")(x)c(+)(x) + c(")(x)B(+)(x) - c^\x)B^{x)}. 

(3.27) 

Thus we realize that the two conditions (3.17) and (3.23) are equivalent under the con- 
straints (3.25). When we define the constrained daughter space V d as a subspace of Vd 
constrained by (3.25), we have 
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Vphys/Vd = Vphys/Vd = U. 

This argument indicates that the condition (3.17) is not alien to QED [19]. 



(3.28) 



3.3 Interpretation of color confinement 

When single-quark states and single-gluon states are not physical, they are unobservable 
and hence confined. Thus the problem of color confinement is settled if we could prove 

Q B \quark) ^ 0, Q B \gluon) ^ 0. (3.29) 

This is a definition of c.c. given in terms of unobservable quantities so that we shall present 
an alternative definition of c.c. in terms of observable quantities in Section 5. 

In order to study the condition on which the relations (3.29) are satisfied we start from 
an identity: 

(A*(x),B b (y)) = (0\T[A-(x),B b (y)] |0) (3 _ 3Q) 
= S ab d fl D F (x -y), 

where a and b denote color indices and Dp is the free massless propagator. This identity 
implies that both A^ and B generate a massless spin zero particle as applied to the vacuum 
state. Thus their asymptotic fields may be expressed as 

A-(x) in = a«(x) + d^x), B b (y) m = (3 b (y), (3.31) 

where denotes the incoming gluon field and x an d P are the incoming fields of the 
massless spin zero particle. Here, we have assumed the validity of the Lehmann-Symanzik- 
Zimmermann (LSZ) asymptotic condition [20] or its suitable modification so that we can 
relate a field operator to a particle state through its asymptotic field. Since QCD is infested 
with infrared singularities an infrared cut-off is introduced to validate the asymptotic 
condition, and only after confinement of colored particles it can be lifted safely for the 
system of hadrons. This is related to the Meissner-like effect discussed in Subsection 6.3. 
The asymptotic fields introduced above satisfy the relations, 

8^ = 0, D 2 a«=0, (3.32) 

n x a = a ■ (3 a , nf3 a = 0, (3.33) 

and 

( X a (x),P b (y)) = -5 ab D F (x-y), (3.34) 

(P a (x),(3 b (y)) =0, (3.35) 
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Now we shall study the implications of c.c. in the properties of the asymptotic fields. 
Assume that is BRS invariant, namely, 

5al = 0, (3.36) 

then 

Q B \gluon) = Q B a£|0) = (-»)<5a£|0) = 0, 
so that a single-gluon state is physical and hence observable. Therefore, c.c. implies 

6a$ + 0. (3.37) 

Now 

= d^ in + g(A li x c) a ' in . 1 ' ; 

When the second bilinear term on the r.h.s. of Eq. (3.38) is absent, there is no asymptotic 
field of unit spin in this expression and we are led to Eq. (3.36). This is really the case in 
perturbation theory. Therefore, a necessary condition for gluon confinement is the non- 
vanishing of (A^ x c) m . This means that there must be a bound state between a gluon 
and a FP ghost, and then and only then we have (3.37) and gluons are confined. Similarly, 
quarks are confined when a quark and a FP ghost form a bound state. Thus the problem 
of confinement reduces to that of bound states. 

One of the present authors (KN) and Okada studied the bound state problem by 
making use of the Bethe-Salpeter equations in the ladder approximation and recognized 
that quarks and probably also gluons are confined when two FP ghosts form a bound 
state [21]. The statement that c.c. is a consequence of the formation of the dighost bound 
state was plausible but not conclusive because of the approximate nature of the above 
treatment. This condition was further refined and reappeared later in a simpler form, 
namely, the condition (3.49) to be introduced in the next subsection. It is not difficult to 
show that the formation of dighost bound states is a consequence of this new condition. 
This relationship is clarified in Appendix A. 

3.4 Conditions for color confinement 

By utilizing the conserved current (3.15) we shall introduce a set of Ward-Takahashi 
identities. First, we define the conserved color charge Q a in terms of the color current 
density in Eq. (3.14) by 

Q a = Jd 3 xJg{x). (3.39) 

Then let us consider a colored field <£ a belonging to an irreducible representation R a of 
the Lie algebra su(3) of the color symmetry, 

[$<%), Q a ] = R a a ^{y), , 
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where $ is the adjoint of <& and R a is an hermitian matrix. 

For the quark field R a is given by T a in Eq. (3.5) or A a /2, where A a denotes Gell- 
Mann's matrix. For the color gauge field obeying the adjoint representation we have 
R% c = —ifabc- Then we have a set of Ward-Takahashi identities of the form 



d fl (SSA^(x),^ a (y),^(z)) 



R a ay S\x - y){&(y), &(z)) - R a 50 S 4 (x - z)(<S> a (y), ¥(z)) 



(3.41) 



This identity follows from Maxwell's equations. When represents a color singlet field 
we have R = and the r.h.s. of Eq. (3.41) vanishes identically. 

In order to derive a condition for color confinement we need some preparations. When 
an operator M is the BRS transform of another operator N, namely, 

M = 5N = i[Q B ,N] T , (3.42) 

M is called an exact operator. Thus its matrix element between two physical states \a) 
and vanishes identically because of the definition of the physical states (3.17). 

((3\M\a) = 0. (3.43) 

Therefore, when the expectation value of an exact operator M in a state I7) does not 
vanish, namely, 

( 7 |M| 7 ) + 0, (3.44) 
it is an indication that the state | 7 ) is unphysical, 

Q B \l)^0. (3.45) 

We are going to exploit this fact in order to derive (3.29). For this purpose we shall 
remove in Eq. (3.41). In momentum space denotes the momentum transfer so 
that we differentiate the Fourier transform of Eq. (3.41) with respect to the momentum 
transfer and take the limit of zero momentum transfer. We shall illustrate this procedure 
in QED by starting from the following Ward-Takahashi identity: 

(P " q)^(p,q) = -KSf\p) ~ Sp\q)). (3.46) 

We differentiate this identity with respect to and then take the limit of q — > p and 
obtain 

r >.f) = -^^(p). (3.47) 
This derivation is valid provided that T^{p,q) does not have a pole at (p — q) 2 = 0. 
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Now define the spin zero projection of 55 and denote it by M M , then apply the above 
procedure to a colored field $> a , then we have 

{p,j3\M£\p,a) oc Rp a x (kinematical factor), (3.48) 

where \p, a) denotes a state involving a quantum of <3? Q with four-momentum p. In this 
derivation, however, we have assumed the absence of the massless pole as in QED, and 
this assumption is expressed by [22] 

55 X a = 0. (3.49) 

The derivation of (3.48) in the absence of the asymptotic fields is discussed in Appendix 
B. 

Then we can refer to the argument based on (3.44) and (3.45) and we conclude that 
the quantum of the colored field is confined. For a color singlet field we have R = 
and its quantum is not confined just as in the case of (3.36). Thus we may conclude that 
Eq. (3.49) is a sufficient condition for c.c. in the sense of (3.29). 

At this stage we quote some examples in which this condition is not satisfied. The 
first example is an abelian gauge theory represented by QED. In this case the quanta of 
the FP ghost fields are free so that they cannot form bound states with charged particles, 
and charge confinement cannot be realized. We can also give an alternative explanation: 
In this case we can easily derive 

iS6 X = -0, (3.50) 

and the condition (3.49) is not satisfied. 

Another example is found when a certain gauge symmetry is spontaneously broken. 
In analogy with Eq. (3.46) we may express the identity (3.41) for the quark field in 
momentum space as 

(P ~ q)»V«(p, g) = -igFZp [S^(p, a) - Sp\q, 0)] . (3.51) 

Assume that the gauge symmetry corresponding to the group index a is broken and that 
the quark masses for the colors a and f3 are non-degenerate. In the limit q — > p the r.h.s. 
does not vanish because of the mass difference, 

S^{p, a )-S^{p,(3)^0. (3.52) 

This implies that V£(jp, q) should develop a massless pole. Physically the appearance of 
this pole is a signal that the Nambu-Goldstone boson has emerged in the form 

55 x a + 0, U55 X a = 0, (3.53) 
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which violates the condition (3.49). 

These two cases serve to explain why the electroweak interactions do not confine any 
particles. Indeed, the electroweak gauge theory is constructed on the gauge group SU (2) x 
U(l), but this symmetry is spontaneously broken and reduced to U(l) corresponding to the 
electromagnetic gauge symmetry. Thus we are left with only an abelian gauge symmetry 
and this explains why the electroweak interactions do not lead us to confinement. 

3.5 Super-convergence relation 

Quite independently of (3.49) a condition for gluon confinement has been derived from 
the study of the renormalization constant of the color gauge field. 
Let us introduce the propagator of the color gauge field, 

(A-(x)A b u (y)) = ^abfd'ke^^Dpik)^, (3.54) 

where 

D F (k), u = (V - ^)D(k 2 ) + a^p. (3.55) 
We introduce the Lehmann representation [23] for D(k 2 ) by 

then Lehmann's theorem gives 

Z^ 1 = jdm 2 p(m 2 ), (3.57) 

where Z3 is the renormalization constant of the color gauge field. Oehme and Zimmermann 
[24] have studied the structure of this propagator in the Landau gauge and have proved 
the superconvergence relation 

Z3- 1 = / dm 2 p(m 2 ) = 0, (3.58) 

when the number of quark flavors Nf is less than 10, namely, in the case 1) of Subsec- 
tion 5.3 as we shall see later. This proof is based on the renormalization group (RG) 
method, and we shall come back to this method later. Then it was recognized by one of 
the present authors (KN) and also by Oehme that gluon confinement follows from this 
superconvergence relation [25,26]. 

Thus it turns out to be an important problem to clarify the relationship between these 
two conditions (3.49) and (3.58) and we shall study it in the next section. 
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4 Renormalization Group and Color Confinement 

In the preceding section we have obtained two kinds of conditions for c.c, and in this 
section we elucidate their relationship in order to find the most fundamental condition. 

4.1 Renormalization group 

The infinitesimal operator of the RG is given by the differential operator 

V = p^ + (3(g)^-^ v (g,a)^, (4.1) 

where p denotes the renormalization point, a the gauge parameter, and jy the anomalous 
dimension of the color gauge field. An element of the RG may be expressed as 

R(p) = exp(pV), (4.2) 
where p is the parameter of the RG and we have the composition law of the group, 

R(p)-R(p') = R(p + p). (4.3) 

Let Q be a function of g, a and p, and we define the running Q by 

Q(p) = exp(pV) ■ Q(g,a,p) , . 

= Q(g(p),a(p),p(p)), 

with the initial condition 

Q(0) = Q. (4.5) 

Then introduce Green's function G(pi; g, a, p) and let its anomalous dimension be 7(5, a), 
then we have a RG equation: 

[V + 1 (g,a)]G( Pi ;g,a,p)=0. (4.6) 
Its running version defined by 

G(p) = ex P (pV)-G(jp i ;g,a, f j,). (4.7) 
satisfies the following equation: 

%G{p) = -i(p)G{p). (4.8) 

Its integral is given by 

G(pf,g,a,p) =exp[fgdp'y(p)} ■ G(pi;g(p),a(p),p(p)). (4.9) 
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The RG provides us with the relationship between renormalized and unrenormalized 
expressions. In order to give a finite value to an unrenormalized expression we have 
to introduce a cut-off A, then higher order corrections tend to decrease for momentum- 
transfer beyond the cut-off A. Then we may assume that the running coupling constant 
g(p) tends to the unrenormalized one or the bare one g$ in the limit p — > oo, 

lim g(p) = g 

provided that the cut-off A is kept finite. In the RG approach we usually formulate the 
initial conditions by keeping A finite, then we can take the limit p — ► oo in Eq. (4.9) to 
obtain 

G(pi;g,a,n) = exp[j °° dp'f(p')} ■ G^\pi;g ,a ,oo), (4.10) 
where denotes the unrenormalized Green function and the factor 

Z = eM-Io °dP7(p')]- (4- 11 ) 
gives the renormalization constant of Green's function G. In particular the gluon propa- 
gator D(k 2 ) satisfies 

R(k 2 ; g, a, p) = k 2 D(k 2 ;g, a, p) = 1, for k 2 = p 2 , (4.12) 

since this is precisely the definition of the renormalization point p. Replacing G by R in 
(4.9) we find 

R(k 2 ; g, a, p) = exp[2 $ dp'j v (p')] ■ R(k 2 -g(p),a(p),p(p)). (4.13) 
Then, by putting k 2 = p 2 (p) in (4.13) and referring to (4.12) we find 

R(p 2 (p);g,a,p) = exp[2 /„' dp'^ v (p')}. (4.14) 

Now insert the Lehmann representation (3.56) into the l.h.s. of Eq. (4.14) and take the 
limit p — > oo or p 2 {p) = p 2 exp(2p) — > oo, and we find, for a finite cut-off, the relation 

Jdm 2 p(m 2 ) = exp[2 j^dpVtp')]- (4-15) 
Combing this relation with Lehmann's theorem (3.57), we find 

Z 3 - 1 = exp[2/ oo dpV(p')], (4-16) 
and it satisfies the RG equation 

(V + 2 1V )Z^ = 0. (4.17) 
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In the cut-off theory we first take the limit p — > oo and then A — > oo, but in what follows 
we invert the order of limiting procedures by taking the limit A — > oo first. Thus some of 
the initial conditions introduced in the cut-off theory are not necessarily satisfied as we 
shall see in what follows. 

4.2 Conditions for color confinement and their relationship 

We shall go back to the condition (3.49). Since i5Sx a is a free massless field we may 
assume that its most general form is given by 

i55x a = -C(3 a . (4.18) 

In QED we have C = 1 and confinement requires C = 0. Now we have to study how to 
determine the coefficient C, and we start from the following relation based on Eq. (3.34): 

(i55 X a (x), X b (y)) = C5 ab D F (x - y). (4.19) 

The field x ' ls a complicated asymptotic field, however, and we shall express this relation in 
terms of Heisenberg operators. For this purpose we shall consider the two-point function 

{i66Afa),A$(y)), (4.20) 

then because of Eq. (3.15) the most general form of its Fourier transform can be expressed 
as 

(<V - fe) / ^ 2 p$& + C$±. (4.21) 
The second term corresponds to the contributions of the spin zero massless particles and 
we immediately obtain 

d^i85Al{x),A b v {y)) = i5 ab Cd v 5\x - y). (4.22) 

Furthermore, because of (3.15) the l.h.s. can be cast in the form of an equal-time commu- 
tator, 



S(x - yo)(0\ [i66A%(x), A){y)] |0) = i6 ab Cd j 6 A (x - y), 

U = 1,2,3). 

As has been shown in [27,28] this constant C satisfies an RG equation 



(4.23) 



(V - 2 7FP )C = 0. (4.24) 

In order to study the equal-time commutator (4.23) we shall go back to the unrenormalized 
version, then 

i56A^ = fyA$ + g^{Af ] * ^) + 9o4° ] , (4-25) 
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where = d^Ay — dyA^ and the superscript (0) is attached to unrenormalized expres- 
sions. When we insert this expression into (4.23) we find that the only surviving term is 
given by 



S(x 
S(x 



yo)(o\ 
yo)(o\ 



i55Al{x)^\A){y)^ 
^ (x)(°),4(y)(°) 

-y), 



|o) 



(4.26) 



where is a parameter which depends on the normalization of the gauge field. In the 
unrenormalized version 
the unrenormalized version 



1 and in the renormalized version a 



7-1 



So we find in 



C(°) = a(°). 

Therefore, in the cut-off theory the boundary condition for C is given by 



(4.27) 



lim(C(p)-a(p)) = 0. 



(4.28) 



The RG equations for C and a are given, respectively, by 



(4.29) 



J^a(p) = -2j v (p)a(p). 



(4.30) 



7i?P denotes the anomalous dimension of the FP ghost fields. The formal solution of the 
RG equation for C{p) satisfying the boundary condition (4.28) is given in terms of a(p) 
by 



C(p) 



a {p) ~ 2 L°° dp \lv{p) + ~tFp{p'))a{p) 



x exp 



-2/; dp"- 1FP { P " 



(4.31) 



After obtaining this formal solution we then let the cut-off A tend to oo. As a consequence 
the boundary condition (4.28) is not necessarily satisfied since we are changing the order 
of limiting procedures. 

Assume that the super convergence relation (3.58) is satisfied so that a(p) vanishes 
identically, then we conclude from Eq. (4.31) that C = C(0) vanishes also. In other 
words, the condition (3.49) or C = follows from (3.58). Therefore, the superconvergence 
relation Z^ 1 = is the most fundamental condition for c.c, and it will be referred to as 
the CCC. This condition was first recognized as the condition for gluon confinement in 
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the Landau gauge [27] since Z 3 had been known only in this gauge, but later this result 
was extended to other gauges and also to other colored particles [28]. 

It should be mentioned that although we have started from (4.28) the difference C(p) — 
a(p) does not necessarily vanish since C{p) and a(p) satisfy different RG equations. In 
particular, C — Z^ 1 = C(0) — a(0) represents the contributions of the so-called Goto- 
Imamura-Schwinger term [29,30]. When Z^ 1 = 0, however, this term should also vanish 
as we have remarked above [31]. We also have 



d(x Q -y )(0\ [A-(x),Ai(y)\ |0> 
-iS ab 5 fll/ Z^ 1 5 A (x - y) = 0, 



(4.32) 



and hence the vanishing of the Goto-Imamura-Schwinger term seems plausible as has been 
proved otherwise. 

We conclude that c.c. is realized when Z^ 1 = 0, although originally only gluons were 
considered to be confined. This CCC implies C = 0, but we already know that C ^ for 
abelian and broken non-abelian gauge symmetries. Therefore, in these cases Z^ 1 cannot 
vanish either. 



5 Realization of Color Confinement 

In the preceding section we have obtained a simple but generic condition for c.c, then a 
natural question is raised of how to evaluate the renormalization constant Z3 so that we 
know precisely in which theories c.c. is realized. It so happens that Z% can be evaluated 
exactly in QCD and we can readily check the CCC [28,32], but before entering this subject 
we have to discuss a more fundamental subject. 



5.1 Gauge-independence of the concept of color confinement 

In QED the renormalization constant Z3 does not depend on the gauge parameter and is 
hence gauge-independent. It is not the case in QCD, however, and it is a function of the 
gauge parameter a and the gauge coupling constant g. Then, what is the significance of 
the CCC since Z^ 1 might vanish in certain gauges but not in others? 

In order to examine this question we first stress that the concept of color confine- 
ment is gauge-independent. When the condition (3.17) is satisfied, the only observable 
particles are hadrons, namely, color singlet bound states of quarks and gluons and they 
are represented by color singlet, and hence BRS invariant, composite operators when the 
LSZ reduction formula is applied. Then the only observable S matrix elements are the 
transition amplitudes among hadronic states. 

Bearing this in mind we introduce the concept of the equivalence class of gauges [28,32]. 
Let us consider a class of Lagrangian densities {C a } representing a gauge theory such as 
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QCD. Assume that all the members of this set are BRS invariant, 

5C a = 0, (5.1) 

and further that the difference between any two members are exact, namely, 

AC = C u - d = SM, (5.2) 

then this set {C a } is called an equivalence class of gauges. Lagrangian densities corre- 
sponding to different choices of a in (3.1) belong to the same equivalence class. 

We introduce Green's functions in two gauges of the same class, then they are related 
to one another through the Gell-Mann-Low relation [33]: 

(A( Xl )B(x 2 ) • • .) n = (A( Xl )B(x 2 ) • • • ea^iAS))/, (5.3) 

where A(xi), B(x2), ■ ■ ■ are local operators, and 

AS = J d^xAC = J d 4 x5M. (5.4) 

In Eq. (5.3) we assume the convergence of the series expansion in powers of AS. 
In particular, when all the local operators are BRS invariant, 

8A = SB = ■ ■ ■ = 0, (5.5) 
we have, as a result of Eq. (3.43), the equality 

(A( Xl )B(x 2 ) ■ ■ -) n = (A(x 1 )B(x 2 ) ■ ■ ■)/ (5.6) 

subject to the convergence of the power series mentioned above. Thus Green's functions 
involving only BRS invariant operators do not depend on the choice of the gauge within the 
same equivalence class. When these composite operators represent hadrons we can apply 
the LSZ reduction formula [20] to such Green's functions to obtain the gauge-independent 
transition amplitudes for hadronic reactions [34,35,36]. 

The unitarity condition for the BRS invariant S matrix between two hadronic states 
| a) and \b) is given by 

(b\a) = (b\&S\a) = J2(b\S f \n)(n\S\a), (5.7) 

n 

and a similar one for SSK Color confinement is realized when the sum over intermediate 
states is saturated by hadronic states alone. Since the hadronic S matrix elements are 
gauge-independent, c.c. in one gauge automatically prevails in other gauges of the same 
equivalence class. This saturation of intermediate states by hadronic states could be 
employed as an alternative interpretation of c.c. since this statement is made in terms 
of observable hadrons alone. We shall come back to this new interpretation later in 
connection with the Meissner-like effect in the present approach. 
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5.2 Evaluation of Z 3 1 

Next we shall proceed to evaluation of Z^ 1 and for this purpose we shall study the RG 
equations for the running parameters g(p), a(p) and fi(p), 

i p g{p)=P{p), (5.8a) 
f p a{p) = -2a(p)j v (p), (5.86) 

f p Ji{p) = p(p). (5.8c) 
First, we shall define their asymptotic values by 

g(oo) = goo, a(oo) = a^, aZ(oo) = oo. (5.9) 
Asymptotic freedom (AF) is characterized by 



goo = o, (5.10) 



and it is realized for Nf < 16. 
By integrating (5.8b) we find 



ln^ = -2 J^dfnv(p), (5.11) 



or 



= exp[2f °° dp^ v (p)] = (5.12) 

Thus evaluation of Z^ 1 reduces to that of aoo. If we identify ctoo with the unrenormalized 
gauge parameter this is a rather trivial relation, but it is not trivial that assumes only 
three possible values 

aoo = —00, 0, a , (5.13) 

where ao is a numerical constant which depends only on the number of quark flavors as 
we shall see in the next subsection. 

5.3 Evaluation of 

Evaluation of has been published elsewhere [17,28,32] so that we shall be brief in what 
follows. 

First, we introduce the (3 function and the anomalous dimension 7^ as series expansions 
in powers of the coupling constant: 
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/%) =5 3 (/?o + /3iff 2 + •••), 



(5.14) 



lv(9, a) = 9 2 {7o(a) + ~fi{a)g 2 + •••), (5.15) 

where 

70 («) = 7oo + 7oia, 

71 (a) = 710 + 7ii« + 7i2« 2 , (5.16) 



The lowest order coefficients are given by 

A) = -^(22-|JV», f . 

Too = -3^(13-1^), 701 = ^>0. { - > 

When fa is negative, namely, when Nf < 16, AF is realized, and we assume it in what 
follows. Then for large values of p we obtain 

g\p)^r P - (6 = -2A,>0) (5.18) 

In order to check the convergence of the integral (5.11) we study the behavior of its 
integrand for large values of p by expanding it in powers of g 2 , and we can easily verify 
the convergence when and only when 

«oo = cxq, (5.19) 

where «o is defined by 

7o(ao) = 7oo + 7oi"o = 0. (5.20) 
When the integral is divergent we obtain 

Qoo = -oo, 0. (5.21) 

Apparently, the CCC (3.58) is satisfied when aioo = — oo because of the sum rule (5.12), 
and we have to check if this case is actually realized in QCD. In order to check which of 
the three possible values of Qqo in (5.13) is realized we have to study Eq. (5.8b) closely 
with the help of AF. In what follows we shall quote only the results. 
Case 1) 700 < 0, a > {N f < 10) 



«oo = < 



a , for a > 0, 

0, for a = 0, (5.22) 
— oo, for a < 0. 
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Case 2) 7o o > 0, a < (10 < N f < 16) 



0, for a > a + h(g 2 ), 
- < «o, /or a = a + /i(g 2 ), 
— oo, for a < ao + h(g 2 ), 



(5.23) 



where 



a = a + % 2 ) = a + S 2 (/*o + h ig 2 + • • •) (5.24) 
is a special solution of the following equation: 

Thus in both cases = — oo is realized when the gauge parameter and the gauge coupling 
constant are properly chosen and consequently the CCC is satisfied. 

A domain in the (a,g 2 ) half-plane corresponding to is denoted by D(a OQ ), then 
this half-plane is covered by D{— oo), D(0), D(ao), where the bar denotes closure. The 
above arguments show that c.c. is realized in D{— oo), but what happens in the other two 
domains? Formally we can apply the argument on the gauge-independence of the concept 
of c.c, namely, for two different values of a we have 

A£ = l(Aa)B-B = -i(Aa)5(c-B), (5.26) 

or 



M = — i(Aa)(c • B), (5.27) 

so that we can utilize the formula (5.6) to show that c.c. prevails in the other two domains. 
In this case Z^ 1 = is a sufficient condition but not a necessary one for c.c. since it is 
not satisfied in the two domains -D(O) and D(olq). The identity (5.3) is valid only when 
the series expansion in powers of AS is convergent. In the present case we have to check 
the convergence of the power series in Aa. Then Eq. (4.9) shows within the framework 
of the RG approach that this implies the convergence of the power series in Aa(p). Let 
us consider, for instance, the case 1), then the line a = is the border between D{— oo) 
and D(ao). Let us introduce a.\ < and 02 > 0, then even when |Aa| = \ai2 — ai\ 1, 
Aa(p) tends to 00 in the limit p — > 00 and the convergence of the power series turns out to 
be doubtful thereby suggesting that a = would be a branch point of Green's functions 
involving BRS variant operators. This, in turn, means that such Green's functions would 
be multi-valued functions of a, Then a question is raised of whether Green's functions 
involving only BRS invariant operators be also multi-valued. When this were the case 
c.c. might be realized only in the domain D(—oo), so that Z^ 1 = would represent a 
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necessary and sufficient condition for c.c, but it has not been clarified yet if it would really 
be the case. 

We may close this section by concluding that c.c. is realized in QCD provided that 
color symmetry is not spontaneously broken and AF is valid. 

6 Consequences of the CCC 

In this section we shall show how the CCC, Z^ 1 = 0, is related to other interpretations 
of confinement based mainly on the dual Meissner effect. Since the starting point of 
other interpretations are quite distinct from that of the present approach it is not easy to 
compare the basic formulations for the purpose of clarifying their relationships. Therefore, 
we shall try to compare the consequences of our approach with those of the others. 

An attempt has been made, however, by one of the present authors (M.C.) and 
Kobayashi [37] to clarify the relationship between the Seiberg formulation [38,39,40] and 
the present super convergence rule in the Landau gauge by studying the criteria for con- 
finement expressed in terms of the (3 function and the anomalous dimension of the gauge 



6.1 An intuitive interpretation of the CCC 

Let us consider a dielectric medium and put a positive test charge inside, then negative 
charges are attracted and positive ones are repelled by it. Therefore, it induces a new 
charge distribution in the medium and the total charge inside a sphere of radius r around 
the test charge is a function of r. It is denoted by e(r) and called the running charge. 

In classical physics the vacuum means the empty mathematical space or the void, but 
in quantum physics the physical vacuum is a kind of medium with a rich structure and has 
to be distinguished from the classical vacuum or the mathematical space. The dielectric 
constant is defined relative to one of the vacua. Now let us regard the physical vacuum 
as a dielectric medium and call the test charge eo = e(0) the bare charge and the total 
charge inside a sphere of a sufficiently large radius e = e(oo) the renormalized charge. 

We then introduce the dielectric constant of the physical vacuum e relative to the 
mathematical one and write down the static Coulomb potential between two electrons in 
two alternative ways, 



We also introduce the renormalization constant Z3 of the electromagnetic field, then we 



field. 




4nr Aner ' 



(6.1) 



find 



,2 _ z^p 2 
— ^3 e > 



(6.2) 
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Thus we are led to 




(6.3) 



At the end of Section 2 we have discussed the two-phase structure or the emergence of the 
Meissner effect in the limit e — > 0, and now through Eq. (6.3) we find that this condition is 
equivalent to the CCC. Furthermore, we find that this condition along with (6.2) implies 
eo = or AF. In the extreme case of e = 0, a small test charge would attract an unlimited 
amount of like charges around it leading the system into a catastrophic state of infinite 
charge. Nature would take safety measures to prevent such a state from emerging, and a 
possible resolution to avoid it would be to bring another test particle of opposite charge. 
The total charge of the whole system is then equal to zero and charge confinement would 
be realized. 

The dielectric constant of the vacuum e is larger than unity, however, as a consequence 
of the screening effect due to vacuum polarization. Quantum mechanically Z^ 1 is larger 
than unity because of the contributions of charged particles of positive-definite metric so 
that the condition Z^ 1 = could be realized only in the presence of charged particles of 
indefinite metric. They do not appear in QED unless they are put by hand just as the 
magnetic monopoles in Section 2, however, so that electric charges cannot be confined in 
QED. Indeed, this is a conclusion repeatedly drawn, but QED is a good laboratory for 
Gadankenexperiment, however fictitious, to illustrate the mechanism of confinement. 

6.2 The CCC and the linear potential 

In this subsection we shall give an intuitive argument on a possible connection between 
the CCC and Wilson's area law in the lattice gauge theory [6]. 

Wilson has formulated the criterion for confinement in terms of a loop correlation 
function defined by 



where P stands for the path ordering. W creates a tube of flux, along the path C, of 
strength equal to one quark color charge, a quark-antiquark pair is attached to both ends 
of the open path or the path is closed. The average value of this expression is defined in 
terms of path integrals in Euclidean space, and for a large closed path C, its average value 
is asymptotically proportional to 



W[C\ = TrPexp(igf c A li dx li ), 



(6.4) 



(W) oc exp(— perimeter) 



(6.5) 



or 



(W) oc exp(— area). 



(6.6) 
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The perimeter and the area denote those of the path C, respectively. Wilson's criterion 
for confinement is the realization of the latter case called the area law. 

In evaluating the loop correlation function (6.4) we shall choose a rectangular contour 
with a temporal extension T and a spatial extension R, and the effective potential between 
a heavy quark and a heavy antiquark is given by 



V[R] = — lim — \nW[R x T\. (6.7) 
When the area law (6.6) is obeyed we have the following linear potential at large distances: 



V[R] = aR. (6.8) 

The linear potential is no longer valid, however, when a quark-antiquark pair can be 
created from the vacuum, since it is energetically more favorable to split the tube of flux 
between the heavy quark pair thereby attaching a light quark and a light antiquark pair 
to the two split ends than to stretch the string indefinitely. We shall realize this point 
later in this subsection. 

In evaluating the expression (6.7) we shall introduce a cluster expansion as given by 

In (TrP exp(ig § c A^dx^)) 
= -b 2 §c§c(M x )My))dx li dy„ + --- • 1 ' ; 

Since we are taking the limit T — > oo for the path C = R x T we may keep only the 
temporal path, and we obtain 



V[R] <x f dx D F (x) 00 , (6.10) 

where Dp{x)^ u is the gluon propagator in Eq. (3.54), and we may insert the representa- 
tions (3.55) and (3.56) into the integrand. 

When the CCC is satisfied we may express D(k 2 ) as 

D(h?)=Sdm \^l u y , (6.11) 

where 

T {m 2 )= jf dM 2 p{M 2 ), (6.12) 

with 



r(0) = r(oo) = 0. (6.13) 

Now we switch to a rather crude approximation since the concept of the static potential 
itself is phenomenological in nature. Actually r = is a singularity in the potential. 
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Although the CCC suppresses the 1/r singularity we are still left with a mild singularity 
at r = 0, and this forces us to introduce a crude approximation. Let us assume that |r(m 2 )| 
has a maximum at m 2 = k 2 , then because of asymptotic freedom gauge interactions are 
stronger at lower energies than at higher energies and k represents a typical mass of 
unconfined systems sharing the same set of quantum numbers with a single gluon. Thus 
we may further assume that k is of the order of or even smaller than the pion mass. Then 
we may approximate (6.11) by a smeared dipole type propagator 

D(k 2 ) « {k .J_ ie) > . (6.14) 

Then the static potential V qq -(r) between a heavy quark and a heavy antiquark can be 
evaluated by utilizing (6.10) as 

V q - q {r) ex l e — 

^-r,forr<^ 1/k. { } 

Thus we have a linear potential in a limited range. The first constant term may be included 
in the self energies of the heavy quarks [41,42]. An unlimited extension of a linear potential 
to large distances would lead to too strong a van der Waals force between hadrons [43] 
so that the potential should cease to be linear and fall off at large distances as remarked 
in the beginning of this subsection. In this connection it should be mentioned that the 
approximate linear potential is a result of but not the cause of color confinement as we 
have seen already. 

The qq interactions discussed here are generated by exchanging colored objects and 
has a relatively long range. This indicates, however, that vortices connecting a quark- 
antiquark pair should have a finite length of the order of or less than 1/k. As we shall see 
in the next subsection the hadron interactions are generated by exchanging color singlet 
hadrons and are of the short-range type. 



6.3 Hadron interactions and the Meissner-like effect 

In the Subsection 6.1 we have discussed the emergence of the two-phase structure and the 
related dual Meissner effect. This is expected to take place in the limit of the vanishing 
dielectric constant. In real QED, however, it is larger than unity and this scenario fails to 
be realized. For instance, let us consider the interaction between two electrically neutral 
systems. Then, their interaction is given by the van der Waals potential 

V vdW {r)<xr- Q . (6.16) 

This shows that the electric field generated by neutral systems penetrate into the vacuum 
without any sharp cut-off indicating the failure of the dual Meissner effect. 
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In QCD we can elucidate the dynamics of hadrons by making reference to dispersion 
relations. It has been clarified by Oehme [44] that dispersion relations for the scattering 
of hadrons remain applicable provided that confinement excludes quarks and gluons from 
the physical state vector space. This is precisely one of the consequences of the CCC, and 
we can extract necessary information about hadron interactions from dispersion relations 
on the assumption that the complete hadron spectrum can be accounted for by QCD. 

Let us consider the nucleon-nucleon scattering as an example, then the potential be- 
tween them is given by the pole contributions in the crossed channels. The least massive 
hadron that can be exchanged between them is the pion, and the resulting interaction is 
represented by the Yukawa potential, 

Vy(r) oc (6.17) 

where denotes the pion mass. 

When we compare this result with the van der Waals potential we recognize that the 
flux of the color gauge field emerging from color singlet nucleons cannot penetrate into 
the confining vacuum leaving no trace of long range forces and that the penetration depth 
is given by the pion Compton wave length. Thus we notice that the Yukawa mechanism 
of generating the nuclear forces bears a strong resemblance to the Meissner effect in the 
type II superconductor. 

Although we have chosen a rather abstract approach to c.c. on the basis of BRS 
symmetry and asymptotic freedom, it shares essentially the same salient features with 
other approaches in that the vacuum allows penetration of the chromoelectric flux from 
hadrons only by a finite depth into it. 

To conclude we stress that all these characteristic features of confinement are the 
consequences of the CCC, namely, Z^ 1 = 0. 
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Appendix A: Derivation of Dighost Bound States from the 
Condition (3.49) 

In this appendix we shall prove the existence of the dighost bound states on the basis of 
the condition (3.49) by utilizing the properties of the representation of the BRS charge 
[16,17]. 

Let us introduce a complete set of basis {|ej}} in the state vector space V and define 
the matrix 77 by 
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Vij = ( e i\ e j)- (Al) 
Then, due to the properties of the inner product, it is hermitian, 

r/t = v . (A.2) 

Provided that V is non-degenerate we can always choose the set of basis so as to satisfy 
the constraint 

V 2 = 1, (A3) 

which defines the standard form of rj. 

The representation t of a given linear operator T with respect to the given complete 
set of basis is defined by 

T\ei) = Ej \ej)tji, (AA) 
and the hermitian conjugate is defined by 

(k\T\l) = (l\T\k)*, (A5) 

for an arbitrary pair of states \l) and \k). Then the representation of denoted by t 
satisfies the relation 

(r?t)t = rjt. (A6) 

The BRS charge Qb is hermitian and nilpotent as discussed in Section 3, and its 
representation q satisfies 

q 2 = 0, = rjqrj. (A7) 
Now we can introduce the Hodge decomposition of V as 

V = V s + V d + V p . (A8) 
The typical members of these subspaces satisfy the following relations: 

q\s) = 0, q^\s) = 0, for \s) € V s (harmonic), 

q\d) = 0, gt|d) ^ 0, for |d) G V d (exact), (A9) 

q\p) / 0, = 0, for \p) G V p (co-exact). 

In each subspace we introduce a complete set of basis such as {|s«)}, {\di}} and {\pi)} 
satisfying 
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(si\-Sj) — (di\pj) — Sij, , . , 

(di\d j ) = (p i \ Pj ) = 0. {Ji - W) 

With these preliminaries we shall start from an analogue of (3.34), 

(x a (.x),p\y))=5 ab D F (x-y), (All) 
where (3 is the asymptotic field of B. Then, by making use of (3.43) and (3.11b) we find 

(6 X a (x), c b (yY n ) = i5 ab D F (x - y). (A.12) 
When the condition (3.49) is valid we have 

5 X a (x)\0) G V s + V d . (A13) 
However, 5x a (x)\0) has zero norm as seen from 

(0\5x a (x)5x a (x)\0) = (0\S (x a (x)S X a (x)) |0> = 0, (A14) 
where we have used (3.43) for Q B . Hence we conclude 

5 X a (x)\0) eV d . (A15) 

The orthogonality conditions (A. 10) and (A.12) then imply that c b (y) m |0) should involve 
a non- vanishing co-exact component and hence 

5c b (yY n \0) = -\g (c(y) x c(y)) b ' m |0) + 0. (A16) 
This establishes the existence of the dighost bound states. 

Appendix B: Proof of the Condition (3.45) for Colored Par- 
ticles 

When the condition (3.49) is satisfied we can prove (3.44) and hence (3.45) for colored 
particles. For this purpose we have to derive (3.48) by applying the reduction formula to 
the identity (3.41). In the absence of the asymptotic fields, however, we have to devise a 
new proof. 

We start from the Lehmann representation [23] for a complex scalar field, 
A' F (z-y) = (0|T U(x)^(y)] |0) 

(B.l) 

= C 2 dM 2 p(M 2 )A F (x-y,M% 
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where m denotes the mass of the quantum of this field, and 

A F (x, M 2 ) = j0j% J ^pe^A^p, M 2 ), (B.2) 

where 

A F {p,M 2 ) = {p 2 + M 2 -ie)- 1 . (B.3) 
There are two alternative possibilities for the nature of the spectral function p(M 2 ): 

case i) 

p(M 2 ) = 5(M 2 - m 2 ) + a(M 2 )6 [M 2 - (m + p) 2 } , p ^ 0. (BA) 
In this case there is a pole term in the propagator on the mass shell. 

case ii) 

There is no pole term due to infrared singularities and p = in this case. 

In the case i) there are asymptotic fields defined by 

<f>(x) = 4> m (x) - J d A yA R {x - y) ■ K y ^(y) 

(B.5) 

= ^ut {x) _ j d 4 yAA ( x _ y) . Ky <t>{y), 

where K y = D y — m 2 denotes the Klein-Gordon operator, and the retarded and advanced 
functions A/j and A^, respectively, satisfy the following equations: 

K x A R (x) = K x A A (x) = -5 4 (x). (B.6) 

In the case ii) the mass shell is a branch-cut in the propagator and there are no discrete or 
normalizable single particle states. Instead of a single particle state, for instance, we have 
a superposition of multi-photon states like a wave packet. Also, as clarified by Bloch and 
Nordsieck [45] , electron scattering is an inclusive reaction since an electron is accompanied 
by an infinite number of infrared photons which escape detection. In this way we have 
only a limited class of observable quantities in the infared singular theories. 

Now we turn to the reduction formula in the case i). The general idea is to replace a 
propagator by a single particle wave function. For instance, 

(0|T [<j>(x)tf(yj\ \0)^(s\^(y)\0). (B.7) 
We shall call this operation 1Z S and it is realized by 
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TZ s (0\T [0(x)0t( y )] |o) = /d 4 x( S |0t(x)|O)(-i)^(O|T [</>(x)tf(yj\ |0> 

= W%)|o>, 

or generally, 

K s (0\T [<p{x)AB ■ ■ ■} |0> = (s\T [AB ■ ■ ■} |0). 



(5.8) 



(B.9) 



When T is replaced by the antichronological symbol T, (—i)K x should be replaced by iK x . 

Amplitudes are not what we can directly observe; so we shall study their absolute 
squares and expectation values. The reduction formula for the former should realize the 
following replacement: 



(0|f [^u)^{v)\ |O)(O|T[0(x)0t( y) ] |o) 

-E s (O|0(n)| S )(,|0t(y)|O) 

- f^ +A)2 dM 2 p(M 2 ) i A+(u - y, M 2 ), 



(5.10) 



where A depends on the experimental condition and corresponds to the Bloch-Nordsieck 
states in the case 2). In the case 1) the reduction formula is given, for A < /x, by 



E s / d 4 xd*y(0\T [. • • \0)iK y {0\<l>(y)\8) 

x{s\tf(x)\0)(-i)K x (0\T[<l>(x)...]\0) 

= E s (0|T[...]| S )( S |T[...]|0). 



(5.11) 



In this case the result does not depend on A and we may take the limit A — > 0. 

By the same token the reduction formula for expectation values is given, typically, by 



/d 4 xA( S |^(x)|O)(-z)^(O|T[^)J(z)0t(y)] |0> x (-i)K x (0\<Kx)\s) 
= (s\J(z)\s). 



(5.12) 



In order to proceed to infrared singular theories we shall discuss QED in what follows. 
We introduce the electron propagator Sf(j>) and put 



(»7 • p + m)S F (p) = a(p 2 ) + ^b(p 2 ^ 



(5.13) 
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Then for the free electron field we have a = 1 and 6 = 0. The behavior of the propagator 
near the mass shell is given by 



a(p 2 ), b(p 2 ) oc (1 + £?y, with c = -(3 - a)£. (£.14) 

This result was obtained by solving Ovsianikov's equation in RG [46]. It shows that the 
mass shell p 2 + m? = is a branch-cut. Now define 



D(d) = — ■ d — m, D{d) = 27 • d — m. 
In the absence of the branch-cut we have 

f d 4 xd*y(s$(x)\0)(-i)D(d x )(0\T [p(x)J{z)^(y)\ |0> 

x(-i)D(d y )(0\m\s) 
= (s\J{z)\s). 



(B.15) 



(5.16) 



In special cases we know the expectation values (J) a priori. For instance, the electric 
charge of the wave-packet-like state \s) is known to be equal to — e, and J should be 
replaced by the space integral of the time component of J M , 



Q = Jd 3 zJ (z), 



(B.17) 



where represents the 4-dimensional electric current. In this case (Q) is independent 
of the parameter A indicating the cancellation of the branch-cut between the propagator 
and the vertex function provided that they are related to one another through the Ward 
identity (3.47). 

Since we have two electron propagators in the Feynman diagram corresponding to 
(tpJ^tp), we have to eliminate one of them by multiplying T>(d) instead of D, where 



V{ip)S F {p) = -1. 

By taking the limit A — ► we have 

f d*xd*y(s®(x)\0)(-i)V(d x )(0\T [p{x)J^z)^{y)] |0) 

x(-i)D(d y )(0\^(y)\s) 
= (s\Mz)\s)C. 



(5.18) 



(5.19) 
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Since the normalization of the wave-packet-like state does not guarantee the proper nor- 
malization of the " on-shell- limit state", we have introduced the constant C. It has to be 
determined by utilizing the known value of (Q). In QED we have 

C=l + (3-«)^. (5.20) 

In the absence of the asymptotic fields we can apply the same procedure to QCD and 
obtain the result (3.48) from the identity (3.14) provided that M M represents the spin zero 
projection of 55A^ [27,28] and (7) the on-shell-limit of a single quark or a single gluon 
state in the sense of Bloch and Nordsieck. 
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